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ABSTRACT  Starting  out  from  the  basic  principles  of  quantum 
mechanics,  the  dynamics  equations  for  stimulated  Raman  scattering 
(SRS)  under  semi-classical  theoretical  conditions  are  obtained. 
Theoretical  verification  is  done  of  the  fact  that  classical 
descriptions  are  only  appropriate  provided  the  distance  between 
upper  and  lower  energy  levels  is  relatively  great.  Numerical 
simulations  are  done  of  several  types  of  effects  associated  with 
stimulated  rotational  Raman  scattering  (SRRS)  when  strong  lasers 
pass  through  the  atmosphere.  In  these  are  included  influences 
associated  with  Raman  rotational  properties,  SRRS  threshold 
values,  and  molecular  "memory"  on  pulse  series  transmission. 

KEY  WORDS  stimulated  rotational  Raman  scattering  Semi- 
classical  theory  Transforming  threshold  Pulse  series 
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I  INTRODUCTION 


In  research  on  stimulated  Raman  scattering,  particle  systems 
are  normally  described  by  the  use  of  classical  resonance  sub¬ 
equations  [  1  ] 


I2L 

dt 


+  re’= 


la 


4cojM 


ElE, 


(1) 


As  far  as  stimulated  vibratory  Raman  scattering  (SVRS), 
which  occurs  in  molecules  under  outside  field  effects,  is 
concerned,  opting  for  the  use  of  classical  descriptions  discussed 
above  can  be  understood.  It  has  already  been  clearly 
demonstrated  that,  when  strong  laser  light  passes  through  the 
atmopshere,  SRRS  is  primarily  involved.  With  respect  to  this 
type  of  molecular  transition  between  rotation  states,  there  is  no 
corresponding  classical  model.  There  are  also  problems  awaiting 
investigation  as  to  whether  or  not  it  is  possible  to  make  use  of 
classical  resonance  sub-equation  (1)  to  apply  descriptions. 

From  the  fundamental  principles  of  quantum  mechanics  dp /dt=  -i[H  ,p]IJi 
and  Tripp)  •  this  article  sets  out  to  solve 

particle  field  Bloch  equations  and  field  polarization  p. 

Moreover,  fields  themselves  still  obey  classical  wave  equations. 


2  DYNAMICS  EQUATION  SETS  ASSOCIATED  WITH  SEMI -CLASS I CAL  THEORIES 
OF 

STIMULATED  RAMAN  SCATTERING 
2 . 1  Particle  System  Bloch  Equations 

Assume  that  particle  systems  possess  the  three  energy  levels 
shown  in  Fig.l.  As  far  as  lasers  in  TEMoo  are  concerned. 


I 


n. 


|3> 


Et,  =  eLCOS\j/i_-  (e*‘-+ 
iff  1=^(0  + 


(2) 


Fig.l 
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|2> 

A, 


|1> 


Under  these  effects,  Raman  scattering  will  occur.  The  scattered 
light  is 
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£s  =  esCOSiAs=  (e‘*’+e-‘*‘) 

(Illegible)  system  Hamilton  quantities  are 

^  =  ^0+^, 


Tico, 

0 

0  - 

-0 

0 

^0  = 

0 

Tico, 

0 

H,= 

0 

0 

I '5 

-  0 

0 

Ticoj  - 

n 

0  - 

(3) 
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Dipole  transition  ^  a  /^a^a  ’  t  )»  ^  Dipole 

transition  matrix  elements  /ii.= -<  l|qr|3>,pj= -<2|^r]3>. 

The  atmosphere  is  primarily  composed  o±  dual  homonuclear 
molecules  such  as  N2,  02,  H2,  and  so  on.  Dual  homonuclear 
molecules,  in  and  of  themselves,  have  no  permanent  dipole 
moments.  Transitions  associated  with  molecular  vibrating 
electrons  still  maintaining  base  states  invariablly  do  not 
produce  induction  dipole  moments  either.  As  far  as  the 
production  of  dipole  moments  is  concerned,  it  is  necessary  that 
they  be  accompanied  by  stimulated  electron  states — under  the 
effects  of  electric  field  incidence,  electron  charge 
distributions  give  rise  to  relative  shifts.  This  is  nothing  else 
than  to  say  that  | 3>  states  belong  to  the  class  of  excited 
electron  states.  With  regard  to  SRRS,  | 1>  states  and  |2>  states 
should  belong  to  two  rotation  states  in  the  same  vibratory  state. 

By  solving  equation  dp/dt=  p= -i[H  ,p]/Ti 
satisfying  density  calculations  p=\il/><[l/\ 


Pu=  {e‘*‘-+e  '*•-)- d'^^Pu  («*'' 

Pii  Pii  *  [  WsP23  d^Piz )  (c  e  ) 

-  £i.idi.Pn-  d[p,y)  €-•*•■)]  (5) 

P33~  *[  (d[_Ei^Pi^~  d^Eipf^)  (e'*''+ e  '*'■) 

+  {ds£sP23-  ds^sPu)  ®  ) 


In  this,  d^=  ii:,/2Ti ,  (a  =  L,5),  ^  Use  6i  to  stand 

for  normal  6i>  phases.  Moreover,  note  6L  =  63  -  51  =  513  ,  6S  = 
53  “  62  =  623,  and  51ambda  =  62  —  51  =  612. 

From  similar  equations  associated  with  p'u  • /'23 
,  solutions  are  done  for  | 13  and  |23  (approximate  solutions). 
Again,  substituting  into  equations  (4),  (5),  and  (6),  high 
frequency  terms  are  eliminated,  and  one  obtains  [2] 


P22 -Pu  =  >■[  d^A£‘*^+ ^23 )  -  d^^L -  cr.j )] 
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(8) 


^3=0 


Here  are  introduced  the  signs  d  =d^e*‘-  {a  =  L  ,S  ,X),  £2^-n2i“  <^i2’ 


Q,  0.  .—  {d\z'i  d}E})IAi 


p^Picos/Ai 


called  drift  frequencies) 


)=1+ - 1 -  ,il,,=  {(o^-cOs)t  -  (A-,-^j).-+(<l>L-<Ds)  =  C02r-Afcz  +  A<D  . 

l  +  (2a)s/Aj.)  ^  .  X  n 

Let  •4'  =  ff|2+ 0^12  >  <^12  '  ’  ^3~ 

P22-P11  ;  ,  co.=  ip-2)a,sin\l/i  ,  (Oj=Cli  + 

Using  (Rl,  R2,  R3)  to  be  the  three  components  constituting  vector 
R,  and  using  (Rl,  R2,  R3)  to  be  the  three  components  constituting 
vector  CO,  ,  one  then  has 


i?=  <o>(  R 


(10) 


That  is,  during  stimulated  Raman  scattering,  particle  systems 
satisfy  Bloch  equations.  Normally,  Qe  «  filambda  .  R  is 
primarily  rotational  around  3  axses. 

Doing  transforms 

rcos^/^j  — sini/^2 
Rj  =  sinij/i  cosrf/^  0 

-rJ  L  0  0  1-1  Lw^- 


/195  Leaving  out  high  frequency  terms,  one  obtains 


(11)  ,  . 

The  three  components  of  vector  R;  are 
.  The  three  components  of  the  vector  2.5  are  (-Q, ,  0 

The  rate  of  time  change  associated  with  Rlambda  will  be  ' 
much  slower  than  R. 

This  article  adopts  P  =  1.  At  this  time,  it  is 
equivalent  to  adopting  E  =  —  e  e‘*’ ,  ia=L  ,S).  ’  This  is 

nothing  else  *  2  *  ’ 

than  saying  that  rotational  wave  approximations  in 
resonance  optics  are  also  appropriate  for  use  with  regard  to  SRRS 
processes . 

Taking  stimulated  particle  attenuations  and 
phenomenologically  introducing  them  into  motion  equations,  one 
has 


«2=  -  (A,-A4>)r,-ru2 


(12) 
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5;=  (A,- A<1>  )uj  +  Q,w^— 


(13) 


(14) 


In  these,  1/ri  =  tl  is  longitudinal  relaxation  times.  1/r  =  t2 
is  transverse  relaxation  times,  we  is  an  equilibrium  value  for 
particle  number  reversal  when  there  is  no  field.  We  =  Qsf  +  At. 
At  is  Fig.l's  detuning  amount  [3]. 

2.2  Field  Equations 

Electric  field  Ea  (a  =  L,S)  satisfies  wave  equation 

,J1  __L  JL  = 

^  dt^  *  cj 

No  is  particle  number  density,  p  =  Spa  is  a  particle 
polarization.  No  p  is  macropolarization  strength,  c  =  co/n  is 
the  speed  of  light  in  media,  n  is  index  of  refraction  (Here, 
chromatic  dispersion  is  ignored) .  co  is  the  speed  of  light  in  a 
vaccuum. 

According  to  the  basic  principles  of  quantum  mechanics,  5 
p=Tr{pfi)  .  Solving,  one  gets 

Ps  Ps  *^23  ) 


In  this,  -  =  =  ,  /  that  is,  . 

Taking  polarization  and  dividing  it  into  two  parts — isophase 
and  nonisophase 


(16) 


Let 

i?,i  =  o-.  +  o-;  ,  R^=  -i  (cr,-o-;). 

.  In  this 

(  > 

Sa  =  L  H 

s  ^23  » 

'^a=sm 

Between 

and  (»,  ,13, 

) 

the  transforms  below 

are  satisfied 

“ -R,.!  r  cosiA,  -sin^/r  n  ru  n 
_  i?^jJ  L  siniA^  J  L  u,  J 

Making  use  of  al3  and  a23  as  already  solved  for,  one  obtains 
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(17) 

.  [dgS^Ux  Pn)] 
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(18) 

1  7 

(19) 

(p33  ~  P22 )] 

(20) 

1  r 

^5=  ^  d^e^v^ 

Al 

Particle  number  conservation  requires  p33  +  p22  +  pll  =  1. 
From  equation  (9)^  one  obtains  p33  =  Cl  .  Therefore,  p22  +  pll 
=C2  .  Cl  and  C2  are  constants.  It  is  possible  to  see  that, 

during  ,  ,  .  . 

stimulated  Raman  scattering  processes,  transitions  only  arise 
between  state  pl>  and  state  p2>.  The  arrangement  of  state  p3> 
particles  is  maintained  invariable  throughout.  p3>  is  designated 
as  intermediate  state.  Cl  and  C2  are  determined  by  initial 

It  has  been  pointed  out  before  that,  with  regard  to  dipole 
transitions,  p3>  is  an  electron  excitation  state.  The  energy 
difference  between  it  and  a  base  state  is  approximately  lOeV. 
Therefore,  Cl  =  O.  C2  =  1.  Because  of  this,  one  solves  to  get 


Pi3-Pl2=  -  ('Va+D 


P33-Pll  =  y  ('Vi-1  ) 


(21) 

(22) 


Taking  equations  (17)-(22)  and  substituting  into  equation  (16), 
one  gets 

Pl=Pl{  (23) 


Pl{  ("'!.■'■  1  )]cosV'i+ ^  ^ 

Okl  i^L 


=-  (a).+  <i>J'p. 


(25) 


Taking  equations  (23)-(25)  and  substituting  into  equation  (15), 
comparisons  of  cos^,  and  coefficients  arrive  at 

equations  associated  with  slow  change  amplitudes  and  slow 

change  phases  of  fields  (consider  a)^»<i)^, 


) 


dz 


d  1  d 

Oz  C  ot  ^  ^  ^  ^ 


(26) 
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(27) 


dt 


)<I>L=-Pd 


»x- 


Ms 


(w,-  1  )] 


dt 


) P^E^V  ^ 


(w,+  l)] 


(28) 

(29) 


In  this. 


P.= 


clTiLi^ 


ia  =  L  ,S) 


3  CLASSICAL  STIMULATED  RAMAN  SCATTERING  DESCRIPTIONS  AS  WELL  AS 
THE  APPROPRIATE  CONDITIONS  FOR  USE 

Let  ^  £,=  £.e-*-  (30J 
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g  =  (m^+ ii;^)e'*®  (AO  =  <1)^^— (tj)  (31) 

If  one  adopts  h-.  =  w,  (w^  =  0 ) ,  ,  then  the  fixed 

solution  equation  set  becomes 

+  7  ^  )i,=  -iPd£sQ--^  (»v.-i)£j 

^17  ■'■■7  "7  -iPsl^LQ'  +  -^  (w.+  l)£s] 

[  4-  +  (r  +  .A,)]  Q'=-  ik,w,E^£, 

k,=  AK.  ^ 

In  this,  2TiA^  .  Normally,  one  has  Ae  «  r. 

we  =  -1,  then  equations  (33)  and  (34)  change  to  become 

^  "7  ^  .  (k,  =  fi,) 

(-7  +r)e’  =ik,ElE, 

When  pump  light  does  not  attenuate,  equations  (35)  and  (36) 

form  a  closed  equation  set.  It  is  nothing  else  than  molecular 
Raman  scattering  (Stokes  wave  and  photon  phonon  wave)  coupling 
equations.  Q  is  simple  positive  molecular  vibration  coordinates, 
also  designated  as  photon  phonon  waves.  In  reality,  from 
equation  (31),  it  is  possible  to  see  that  Q  reflects  changes 
associated  with  particle  systems  (“i-i'j)  •  At  the  same 

time,  it  also  reflects  changes  in  field  phase  (|^).  So  long  as 
field  phase  differences  |4>  are  maintained  constant,  Q  then,  and 


(32) 

(33) 

(34) 

If 

(35) 

(36) 


8 


only  then,  belongs  purely  to  particle  “coordinates". 

Key  conditions  for  transitioning  to  classical  descriptions 
are  w^=w,=  -l,  ,  that  is,  particle  lay  out  numbers 

associated  with  transition  upper  energy  levels  are  zero  (p22  = 

0) .  The  ratio  associated  with  first  molecular  vibration 
excitation  states  and  base  state  lay  outs  are  [4]  p22  /  pll  = 
1.4x10-5.  Therefore,  classical  descriptions  are  capable  of 
giving  good  results  with  regard  to  molecular  vibration  Raman 
scattering.  With  regard  to  arrangement  differences  associated 
with  two  rotational  energy  levels  on  vibration  base  states,  they 
are  [5]  we  =  -0.035.  It  is  possible  to  see  that  classical 
equations  will  have  relatively  large  errors  with  regard  to 
describing  SRRS.  However,  provided  corrections  are  made  to 
|N=No|we|,  classical  equations  are  still  capable  of  giving  good 
results . 


4  CALCULATION  RESULTS 

Fixed  solution  equations  are  simultaneous  equation  sets 
composed  of  7  equations — (12)— (14)  and  (26)— (29). 

4 . 1  Parameter  Determination 

When  material  reaches  a  stable  state  w^=0),  / 

it  is  possible  to  solve  and  obtain 

„  = - - ! -  (37) 

r^+a'/r,)a]+iA-A<t>y 

(Illegible)  the  equations  above  and  substituting  into  equation 
(28),  it  is  possible  to  obtain 

d  1  5  , 

^  dz  c  dt 


In  this  /  = ,  (a=L  ,S); 

4n 

_ Sn-NA-wJr _  y 

criAs[r'+ (r/r,)£2j+ (A,-A*<&)^]  ~  m,. 


(Illegible)  stable  state  gain  coefficient.  The  setting  up  of  the 
second  equality  in  the  equations  above  is  because  Qe  «  r  and 
(A,-A*<D)  «  r.  /198 

Kererence  ( 5 )  gives 

647rX(-w.)  (/■+l)(/+2)  , 

^  157iMsAv,c^  (2;+ 1 )  (2y-l- 3  ) 


Therefore, 


,  ^2^  _L  (/•+  1)0+2)  , 

'  M,  15  (2y+l)0-+3) 
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;r  =  r,  ,  Ar=0  ,  A^=  10- 

(photon  phonon  wave  number )  .  Adopting  a's  .  -  cm  ,  a,=  A, -  ,  vl^  ^ 

is  photon  phonon  wave  length The  values  are  solved  for  from 
photon  phonon  wave  numbers  V 

Materials  parameters  are  seen  in  Table  1 
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Table  1 


Table  1  The  parameters  of  materials 


N„  (0) 

Voo 

B 

P  0')/10-'-ns"' 

/I0'’cm-' 

/I0-"cm’ 

/cm*' 

6 

8  10  12 

14 

N: 

2.237 

0.693 

2.001 

-0.035 

3560 

3270  3060  2870 

2660 

O: 

1.099 

1.4378 

-0.039 

4.2  Fixed  Solution  Conditions  /-n\_ 

Initial  conditions:  when  t  =  0,  (0) -v,  (0) -0,  (0) -w^  . 

Boundary  conditions:  when  t  =  O,  at  ground  locations  (z  = 

0),  add  a  planar  wave  propogating  along  the  z  direction.  These 
are  Gaussian  distributions  following  along  with  time  changes. 

Wave  lengths  are  |l  .  Half  widths  are  t .  Peak  value  strengths 
are  IL  (0).  Seed  illumination  IS  (0)  is  adopted  as  IS  (0)  =  IL 
(0)  X  10-7. 

4.3  Calculation  Results 

4.3.1  SRRS  Transformation  Characteristics  and  Threshold  Values 
When  strong  lasers  pass  through  the  atmosphere,  there  is  an 
intimate  relationship  between  scattered  light  strengths  IS 
associated  with  transformations  and  initial  incident  light 
strengths  IL  ( 0 ) .  Concrete  calculation  results  are  as  shown  in 
Fig. 2.  When  IL  (0)  >  3.9MW/cm2,  increases  in  IS  following  along 
with  IL  (0)  are  relatively  slow.  When  IL  (0)  =  3.9  MW/cm2  ,  IS 
increases  rapidly  following  along  with  increases  in  IL  (0).  When 
IL  (O)  =  3.9  MW/cm2  ,  IS/IL  (O)  «  1%.  This  pattern  during  strong 
laser  atmospheric  transmission  possesses  universality.  The 


normal  definition 


_ 

h  (0) 


=  v 


Fig.  2 


/i/MW*cin-= 

Fig. 2  Stokes  intensity  Is  as  a  function  of  the 

incident  intensity  (0 )  1 1 


of  IL  (0)  is  SRRS  threshold  value  It.  When  IL  (O)  >  It  ,  emited 
laser  light  IL  will  be  scattered.  Therefore,  during  research  on 
strong  laser  light  passing  through  the  atmosphere,  study  of  It 
for  various  types  of  wave  lengths  |l  and  pulse  widths  t  is 
significant.  /199 

Table  2  gives  threshold  values  solved  for  for  different 
(T, |l). 

Table  2 


Table  2  Threshold  of  different  (t  , 


t/ns 

0.01 

0.01 

0.02 

0.1 

1.0 

2.0 

2.0 

ijtim 

0.35 

1.0 

1.0 

1.0 

l.O 

l.O 

0.4 

/,/MW  .  cm-  = 

322 

928 

464 

98.5 

14.8 

9.8 

3.9 

4.3.2  Pulse  Series  Transmission  Characteristics 

Pulse  series  are  composed  of  multiple  pulses .  The  width  of 
each  pulse  is  t.  The  interval  between  various  pulses  is  t. 
Therefore,  what  is  different  from  single  pulses  about  the 
transmission  of  pulse  series  is  due  to  molecules  possessing 
"memory"  effects.  Speaking  in  concrete  terms,  this  is  nothing 
else  than  a  previous  pulse  later  being  forced  by  stimulated 
molecules  into  exponential  attenuations  as  a  function  of  time 
using  e~^‘  rules.  See  equations  (12)  and  (13).  In  this  way, 

with  t  <  1/r,  when  the  second  pulse  arrives,  the  molecular 
stimulation  which  the  first  pulse  creates  is  still  partially 
retained.  The  second  pulse  is  scattered  by  this  type  of 
stimulated  molecule.  If  t  »  1/r,  obviously,  pulse  series  are  no 
different  from  the  transmission  of  single  pulses. 
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Fig  .3  The  influence  of  “  the  memory  of  the 
molecule'  upon  the  propagation  of 
the  pulse  series,  as  a  function  of  the 
time  begining  with  the  pulse. 


Fig. 3  The  influence  of  “the  memory  of  the  molecule"  upon  the 
propogation  of  the  pulse  series,  as  a  function  of  the  time 
beginning  with  the  pulse.  (1)  Dimensionless  Time 


/200 

What  Fig. 3  graphs  out  is  curves  (  )  associated  with  material 

stimulation  created  by  pulse  series  pulses  with  serial  numbers  1, 
2,  and  10.  t  is  0.02ns.  IL  (0)  =  0.085MW/cm2  .  Fig. 3(a)  is  t  = 
2ns.  This  clearly  shows  that  molecules  have  "memories". 
Stimulations  are  cvimulative.  Fig. 3(b)  is  t  =  20ns.  The  three 
lines  converge  together  presenting  single  pulse  characteristics. 

Due  to  collision  attenuation  coefficient  r,oce-'“  ,  the 
higher  the  atmospheric  level  is,  therefore,  the  stronger  is 
molecular  “memory",  and  the  clearer  accvimulative  effects  are. 

Fig. 3(b)  is  a  case  where  the  altitude  is  10km.  Fig. 3(c)  is  a 
case  where  the  altitude  is  50km  (t=20ns).  Molecules  have 
"memories",  and  stimulations  are  cumulative. 
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5  ANALYSIS  AND  DISCUSSION 


The  driving  term  Qe  associated  with  in  equation  (14) 

IS  nopally  very  small.  Therefore,  w,  deviation 

equilibrium  values  we  are  very  small.  This  is  nothing  else  than 
the  reason  improved  classical  descriptons  [6]  give  good  results. 

Due  to  the  fact  that  the  driving  term  Qe  associated  with  ’ i 
in  equation  (12)  is  far  larger  than  the  driving  term  (A  -  A*<I>)t). . 
associated  with  in  equation  (13),  material  polarization," 

therefore,  is  primarily  characterized  by  nonisophasal  component 


\\\ 

Seen  from  the  view  point  of  equations  (26)  -  (28), 
directly  influences  changes  in  field  amplitudes  '  .  Then 

*1’  and  e,  directly  influence  changes  in  field  phases 
®  .  Field  amplitudes  £.  .  and  field  phases  (t>  do 

not  directly  give  rise  to  coupling.  By  contrast,  coupling  is 
given  rise  to  indirectly  through  materials  systems. 

When  setting  up  computational  plans,  use  is  made  of  this 
type  of  coupling  property  to  make  program  efficiencies  relatively 
high.  In  5  million  operations  a  second  computers,  program  cvcles 
each  only  require  100  /is.  i 
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ABSTRACT  This  article  introduces  measurement  principles 
associated  with  atmospheric  coherence  length  and  isoplanatic 
^^9!®  measuring  instruments  as  well  as  the  basic  structures  of 
instruments .  It  reports  the  results  of  measurements  carried  out 
at  the  Yunnan  Observatory  on  Kunming's  Phoenix  Hill.  The  average 
value  measured  for  ro  is  12.4cm.  Average  value  for  |o  is  10.6 
Mrad.  Compared  to  data  associated  with  temperature  pulsation 

probe  instruments,  the  results  obtained  were  demonstrated  to  be 
reliable.  Finally,  analyses  were  carried  out  with  regard  to  a 
number  of  factors  influencing  measurement  errors  which  can 
possibly  exist  in  instruments. 

KEY  WORDS  Atmospheric  coherence  length  Isoplanatic  angle 
Atmospheric  seeing  Adaptive  optics 
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I  INTRODUCTION 


In  situations  where  use  is  made  of  adaptive  optics  systems / 
it  is  necessary  to  understand  atmospheric  seeing  at  the  places 
and  times  concerned.  The  selection  of  sites  for  observatories 
and  atmospheric  seeing  associated  with  daily  astronomical 
observations  are  also  one  of  the  indispensible  myriad  parameters. 
Quantitative  descriptions  of  atmospheric  seeing  parameters  are 

nothing  else  than  the  so  called  Fried  length  [6]  -  atmospheric 

coherence  length  ro.  This  characterizes  atmospheric  coherence 
properties  along  certain  transverse  distances  in  a  certain 
specially  designated  light  path.  In  the  same  way,  people  are 
concerned  with  coherence  characteristics  between  different  light 
courses  within  certain  ranges  of  distances  or  certain  angles. 

This  then  requires  making  use  of  isoplanatic  zones  or  isoplanatic 
angles  in  order  to  make  descriptions.  In  actuality,  what  are 
made  relatively  more  use  of  are  isoplanatic  angles,  normally 
recorded  as  ©o.  In  the  middle  1980 's,  we  began  to  make  efforts 
in  research  work  relating  to  atmospheric  coherence  [2].  In 
conjunction  with  this,  an  instrumentation  unit  to  measure  ro  and 
©o  has  been  recently  developed  and  manufactured  successfully. 

This  article  is  nothing  else  than  the  carrying  out  of  a  simple 
introduction  to  its  measurement  principles,  instrument  structure, 
and  preliminary  measurement  results.  Finally,  analysis  is  done 
of  error  factors  which  can  possibly  exist. 

2  MEASUREMENT  PRINCIPLES 


2 . 1  Atmospheric  Coherence  Length 

The  definition  of  atmospheric  coherence  length  is 

D  (r)  =  6.88  (r/ro (rad^) 


(1) 


That  is,  when  receiving  optical  system  radius  r  is  equal  to 
coherence  length  ro,  the  value  of  wave  structure  funcion  D  (r)  is 
6.88  rad2  .  At  this  time,  the  properties  of  the  systems  in 
question  reach  optimum  limits.  Due  to  random  propogation  theory 
associated  with  waves  in  media,  in  air-ground  light  courses, 
planar  waves  ro  ,  are  capable  of  being  expressed  as 


ro 


[0A22khtc\l/ 


r  “ 

Cl{h)dh]-^i^ 

J  h' 


(2) 


In  the  equation,  k=2rr/A,.  X  is  wave  length,  tjr  is  the  zenith 
angle.  C\  (h)  is  height  distributions  associated  with  turbulent 
flow  structure  constants,  h'  is  the  altitude  of  measuring 
stations . 

If  it  is  possible  to  measure  C\  (h),  it  is  then  possible 
to  calculate  ro  from  equation  (2).  Barletti  and  others  [3]  used 
this  method  to  obtain  height  distributions  associated  with  ro. 
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The  results  were  the  same  as  r  distribution  models  [4]  we  use/216 


Cl  (/i )  =  Co /i  exp  (-h /Ho) 


Calculation  results  agree.  In  the  equation,  Co  and  Ho  are  two 
set  up  parameters .  Ho  is  analogous  to  the  equivalent  atmospheric 
altitude  at  the  place  and  time  in  question.  Co  corresponds  to 
turbulence  strengths  associated  with  measurement  station 
positions.  However,  normally,  making  use  of  this  tpe  of  method 
in  order  to  precisely  determine  ro  is  very  difficult.  Moreover, 
it  is  also  not  possible  to  obtain  results  in  real  time.  As 
result,  it  is  necessary  to  opt  for  the  use  of  a  type  of  indire 
mSIsur;ment  method.  In  it,  use  is  made  of  the 

wide  SDread  methods — which  make  use  of  arrival  angles  in  order  t 
do  conversions.  The  theoretical  expression  associated  with 
arrival  angle  fluctuation  variences  is 


<7^  =  2.905  sec  i/^ 

a  ' 


Cl{h)dli 


(4) 


a-  n  ic;  i-he  diameter  of  the  receiving  telescope 

it  it  possible  to  obtain 


ro  =  3.18fc 


1/5  ^-6/5 


(5) 


In  earlv  periods,  measurements  associated  with  arrival 
«nales  mSde  Lrof  photographic  time  exposure  methods  in  order  to 
record  star  fluctuations.  From  one  dimensional  directional  time^ 

exDOSures  it  is  possible  to  obtain  angular  variance  C, 

exposures,  luxb  p  equal  to  each  other,  an 

In  order  to  make  ena  ^  4.k=4-  = 

empirical  constant  exists  between  the  two,  that  is. 


(6) 


Ac:  far  as  adopted  b  values  are  concerned,  there  exists  very  great 

Sve  gllduinf become  known  Is  -duel  point  difference  arrivnl 
anale  methods”  or  "difference  imagery  movement  methods  t6] 
?SstSaS!  This  type  of  method  takes  light  coming  from  a  certain 
ralestial  bodv  and,  through  certain  measures,  fonns  two  celest 
SagSs  on  the  receiving  plane.  Respective  solutions  are  made  for 
the^differences  Apci  of  the  two  celestial  images  from  their 
central  position.  Through  the  formula 


=  r  2 /H  0-36  (A/D) '^^-0.242  (X/d)''^]X 
<Apli>  -  <Apci>' 


5/3 


} 


3/5 


(7) 
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ro  is  calculated.  In  the  equation,  f  is  telescope  focal  length, 
d  is  the  central  interval  associated  with  the  two  light  paths  ( 
that  is  also  the  theoretical  interval  associated  with  two  images 
when  there  is  no  fluctuation).  This  type  of  method  is  not 
sensitive  to  such  factors  as  perturbations  in  the  instruments 
themselves,  optical  system  quality,  temperature  effects  on 
telescope  focal  length,  fluctuations  in  star  image  brightness, 
and  so  on.  As  a  result,  it  is  possible  to  hope  to  obtain  very 
good  measurement  precisions. 

2.2  Isoplanatic  Angles 

The  theoretical  expression  associated  with  isoplanatic 
angels  is  - ® 

0o=[2.9O5it'seciA  J  C\{hWUh]-'^'^ 

In  the  same  way  as  ro,  ©o  daily  observations  must  also  go  through 
observations  of  another  physical  quantity  in  order  to  carry  out 
conversions.  The  most  convenient  solution  is  to  use  star  light 
scintillation  methods.  This  is  due  to  the  fact  that  the 
scintillation  method  also  possesses  altitude  weighting 

functions  similar  to  0o  .  Theoretical  calculations  clearly  show 
that,  when  D  =  11cm,  it  is  then  possible  to  make  the  weighting 
functions  associated  with  the  two  equal  to  each  other  [ 7 ] .  In 
conjunction  with  this,  the  relationship  set  out  below  is 
established 

0o=C[log(l  +  <7,V<5>*)]-^'’  (9) 


<S>  is  average  values  associated  with  star  light  signals.  C  is  a 
constant.  Through  large  amounts  of  experimentation,  Eaton  and 
others  [8]  indicate  that  C  =  0.9676.  Equation  (9),  during  weak 
turbulence  and  medium  strength  turbulence,  the  precisions 
possessed  are  very  good.  It  is  only  under  strong  turbulence 
conditions  that  errors  can  reach  approximately  20%.  In  general 
cases,  turbulence  strengths  throughout  the  entire  atmosphere  will 
not  be  in  a  strongly  turbulent  state.  As  a  result,  equation  (9) 
is  normally  effective.  /217 

3  INSTRUMENT  CONSTRUCTION 

In  our  coherence  length  and  isoplanatic  angle  measurement 
instruments,  as  far  as  measurements  associated  with  ro  are 
concerned,  use  is  made  of  dual  point  difference  arrival  angle 
methods.  With  regard  to  measurements  associated  with  0o  ,  use  is 
made  of  star  light  scintillation  methods.  The  schematic  diagram 
is  as  shown  in  Fig.l.  Principal  parameters  are  seen  in  Table  1. 
The  instruments  in  question  break  down  into  two  large  parts — 
optical  receiving  systems  and  data  processing  systems .  The 
former  includes  two  recieving  telescopes  connected  together  with 
completely  uniform  properties,  image  merging  light  paths  and  two 
detectors  (4720  model  CCD  pickup  cameras  and  photoelectric 
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multiplier  tubes).  The  entire  system  is  installed  on  an 
equatorial  type  base.  It  is  possible  to  use  computers  to  control 
step  electrical  equipment  to  make  searches  and  track  in 
declination  and  right  ascension  directions.  The  latter  includes 
video  frequency  cards,  A/D  cards,  and  386  microcomputer  systems. 


Fig.l 


Fig.  1  Schematic  diagram  of  the  r©  and  9q  measurer 


Table  1 


Table  1  The  main  parameters  of  the  r,  and  d, measurer 


Optical  system 

Receiving  aperture 

Z>/cm 

12 

Effective  focal  length 

//m 

5 

Se  para  ton  of  centers 

d/cm 

15 

Visual  angle  of  telescope 

(') 

3 

Visual  angle  of  guiding  telescope 

n 

0,5 

Image  signal 

Minimum  illumination 

E^/U 

0.02 

Signal -to -noise  tatio 

SNR/dB 

52 

Sampling  rate 

R/Hz 

25 

Sampling  time 

t/s 

2.5 

Intensity  signal 

Dark  current 

/o/nA 

50 

1 

Sampling  rate 

R/Hz 

1000 

Sampling  time 

tis 

5 

Mechanical 

Type 

Equatorial  style 

system 

Right  ascension 

Declination 

Tracking  precision 

Weight 

±100  ° 

±45  ° 

I<l200kg 

/218 

The  height  of  the  entire  instrument  is  approximately  1.2m.  The 
length  is  approximately  1.6m.  Right  ascension  and  declination 
coordinates  are  picked  up  from  21  place  photoelectric  coding 
disks.  At  the  same  time  as  numeric  tubes  display  right  ascension 
and  declination  coordinate  values,  they  are  transfered  to 
microcomputer  to  be  recorded  in  order  to  act  as  a  basis  during 
normalized  data  processing. 

CCD  pick  up  cameras  possess  characteristics  of  low 
illumination  and  high  signal  to  noise  ratios.  In  order  to  be 
advantageous  to  the  reception  of  celestial  images  with  relatively 
weak  optical  powers,  the  signals  go  through  video  frequency 
collection  cards,  and  are,  first  of  all,  entered  into  internal 
computer  storage.  The  collection  rate  is  25  frames/s.  Each 
iteration  of  observations  is  2.5  s.  After  that,  solutions  are 
done  for  the  two  central  positions  pli  and  p2i  ,  thereby 
obtainingAp  ci  =  p2i  -  pli  .  Finally,  ro  is  calculated  from 
equation  ( 7 ) . 

When  measuring  0o  ,  computer  processing  inputs  celestial 
image  strength  signals  by  A/D  cards.  Collection  speed  is  IkHz. 
The  sampling  period  for  each  set  of  data  is  5s.  Reckoning  up 
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signal  fluctuation  variance  and  average  value  <S>  ,  by 

using  fomula  (9),  it  is  then  possible  to  calculate  0o  . 

In  order  to  resolve  such  problems  as  inadequate  computer 
internal  storage  capacity  and  excessively  long  calculation  times, 
when  processing  celestial  image  signals,  option  is  made  for  the 
use  of  "window"  methods.  The  size  of  windows  must  pay  attention 
to  all  of  imagery  fluctuation  ranges,  mechanical  adjustment 
precisions,  and  observation  times.  Considering  imagery 
fluctuations,  they  are  generally  not  larger  than  2".  Mechanical 
precisions  are  approximately  4".  As  a  result,  window  height 
should  not  be  less  than  10".  This  corresponds  to  the  64  lines  of 
a  512x512  imagery  card.  Window  widths  should  also  be  at  least 
10".  The  total  of  the  two  imagery  paths  is  20".  That  is,  window 
size  is  64  lines  x  128  columns.  During  observations,  first  of 
all,  adjusting  telescope  right  ascension  and  declination 
coordinates  makes  celestial  imagery  enter  into  the  window  placed 
on  the  top  of  the  screen.  After  it  is  appropriately  placed,  the 
base  then  stops  rotating.  After  the  center  of  gravity  of  images 
distinguished  by  computer  are  placed  at  suitable  locations, 
number  gathering  begins.  Windows  follow  celestial  images  as  they 
shift  downward,  synchronously  tracking.  After  completing  the 
gathering  of  predetermined  amounts  of  data,  calculations  are  then 
carried  out.  From  the  beginning  of  collection  to  the  presenting 
of  results  only  requires  4  seconds.  Utilization  and  operation  of 
measurement  instruments  is  very  simple.  Automation  levels  are 
very  high. 

4  MEASUREMENT  RESULTS 

From  November  to  December,  1992,  we  carried  out  measurements 
to  be  completed  in  40  days  at  the  Yunnan  Observatory  on  Phoenix 
Hill  in  Kunming.  This  observatory  is  located  25‘’N  at  2000m  above 
sea  level .  Due  to  the  fact  that  the  weather  is  sometimes  rainy 
or  cloudy,  only  20  segments  of  useful  data  were  acquired.  The 
time  period  for  each  iteration  of  observation  was  between  18:40  ~ 
24:00.  Through  analysis  of  measurement  results,  it  was 
discovered  that  the  average  value  of  ro  is  12.4cm  (Unitized  to 
the  zenith  location.  The  same  below.).  The  range  of  changes  was 
2-2 8cm.  The  average  value  associated  with  ©o  is  10.6  jurad.  The 
range  of  changes  is  3-27  urad.  Fig. 2  and  Fig. 3  are, 
respectively,  observation  results  for  ro  and  0o  on  the  two  days, 

4  and  5  December.  From  these  Fig.'s,  it  is  possible  to  see  that 
changes  associated  with  ro  are  relatively  smooth.  However,  ©o 
changes  are  still  relatively  large.  This  clearly  shows  that 
changes  associated  with  low  atmosphere  turbulence  are  not  large. 
However,  atmospheric  turbulence  flows  associated  with  relatively 
high  layers  (h  >  3  -  5  km)  are  very  unstable.  Looking  at  this 
overall,  the  atmospheric  seeing  for  Kunming  in  this  period  was 
very  good.  It  can  be  compared  to  excellent  observatory  sites  on 
a  world  scale,  such  as  the  U.S.  Mauna  Kea  and  Chile's  La  Silla. 

During  periods  when  use  was  made  of  coherence  length  and 
isoplanatic  angle  measurement  instruments  to  carry  out 
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observations  of  ro  and  0o  ,  Zeng  Zongyong  and  others  fired 
temperature  fluctuation  space  probe  instruments  at  the  same  time 
to  carry  out  measurements  of  c^„  •  Fig. 4  is  an  actual  example 

measured  at  23:00  on  29  November" 1993 .  In  the  Fig.,  curve  b  is 
the  Huf nagel  model  [ 9 ] .  Curve  c  is  given  by  equation  ( 3 ) .  In 
this,  Co  =  3.5x10-14  .  Ho  =  5400m.  Their  standard  deviations 
relative  to  experimentally  measured  data  are,  respectively,  4.6 
and  0.026.  Coefficients  of  correlation  are,  respectively,  0.70 
and  0.73.  From  this,  it  is  possible  to  know  that  r  models  agree 
relatively  well. 

We  made  use  of  probe  data  and  ro  and  ©o  observation 

data  for  the  same  periods  to  carry  out  comparisions .  In 
accordance  with  equation  (2)  and  equation  (8),  respectively, 
calculations  were  gotten  of  ro  =  14.4cm  and  ®o  =  8.3  /xrad.  These 
two  are  in  basic  agreement.  However,  there  are  some  relatively 
small  differences.  The  causes  for  the  differences  that  exist 
will  be  discussed  in  the  sections  that  follow.  /219 
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Fig.  4  An  e.\ample  of  measured  C\  (23:00,  Nov.  29, 1993) 
a .  measured  data ;  b .  Hufnagel 's  model ;  c.  Eq .  (3 ) 


5  CONCLUSIONS 


The  fact  that  there  exist  differences  associated  with 
measurement  values  and  theoretical  calculation  values  is  also 
seen  in  the  work  of  other  scholars.  For  instance,  Sydney  and 
others  [10]  discovered  ro  observed  values  to  be  4.4cm.  However, 
calculated  values  were  3.35cm.  The  two  differed  by  24%.  The 
differences  associated  with  0o  were  even  larger.  This  is 
entirely  consistent  with  our  results.  In  fact,  it  is  reasonable 
that  certain  differences  should  exist.  The  reason  is  the 
measurement  instruments  themselves.  No  matter  whether  it  is 
measurements  of  ro  and  ©o  or  measurements  of  such  parameters  as 
C;  ,  there  always  exists  a  certain  measurement  error.  As 
random  quantity  statistics,  differences  of  around  20%  are  quite 
small.  At  the  same  time,  theoretically,  it  is  also  possible  for 
certain  flaws  to  exist. 

Speaking  in  terms  of  measurements,  our  imagery  gathering 
speeds  are  not  high  enough.  Sampling  times  are  not  long  enough. 
This  then  creates  frequency  spectra  in  signals  not  being  capable 
of  complete  coverage.  Seen  on  the  basis  of  frequency  spectrum 
fluctuations,  high  frequency  components  (f^lOHz)  are  already 
extremely  weak.  Very  large  errors  are  not  necessarily  created. 
However,  paying  inadequate  attention  to  low  frequency  components 
(f:S0.2Hz)  will  give  rise  to  relatively  large  errors.  As  a  /220 
result,  sampling  times  should  be  increased  later.  Secondarily, 
there  is  only  a  15cm  distance  between  the  centers  of  the  two 
telescopes  associated  with  measuring  instruments. ^  It  may  be 
somewhat  smaller.  On  the  basis  of  Roddier  s  opinions,  when  d/D^^ 
2,  the  accuracy  of  formula  (7)  is  somewhat  higher.  However,  this 
is  not  important.  The  shortening  of  interval  distances  will,  no 
doubt,  influence  low  frequency  components.  However,  on  the  basis 
of  the  Taylor  principles  of  frozen  turbulence  flows,  if  option  is 
made  for  the  use  of  multistage  averaging,  it  will  reduce  in  the 
same  way  the  effects  of  low  frequency  components.  One  possible 
cause  for  the  existence  of  relatively  large  differences  in 
isoplanatic  angles  is  that,  according  to  theory,  D  should  be 
11cm.  However,  what  we  actually  used  was  12cm.  Due  to 
relatively  large  apertures  possessing  relatively  large  aperture 
smoothing  effects,  this  leads  to  signal  fluctuations  becoming 
relatively  weak.  Thus,  relatively  large  isoplanatic  angles  are 
obtained.  Another  possible  cause  which  is  more  important  is 
that,  when  solving  formulae  (2)  and  (8),  consideration  was  not 
given  to  the  influences  of  turbulent  flow  outside  scales.  As  far 
as  problems  in  this  area  are  concerned,  we  have  already  made 
detailed  analyses  [11]  in  another  report.  It  is  not  necessary  to 
repeat  them  here.  We  only  wish  to  point  out  this.  After 
considering  finite  outside  scales,  based  on  a  certain  outside 
scale  height  distribution  model,  it  is  possible  with  regard  to 
the  turbulent  flow  data  discussed  above — to  do  calculations  and 
obtain  ro  =  20.9cm  and  ©o  12  Mrad.  These  numerical  values  and 
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observed  values  are  then  in  quite  good  agreement. 

In  summary,  we  are  able  to  reach  the  following  conclusions. 
Results  obtained  from  coherence  length  and  isoplanatic  angle 
measurement  instruments  are  reliable.  After  making  certain 
improvements,  it  will  then  be  possible  to  obtain  even  higher 
accuracies . 
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ABST^CT  This  article  introduces  a  type  of  laser 
conditioning  technology  associated  with  raising  high  power  laser 
reflector  mirror  coating  damage  threshold  values.  In  conjunction 
with  this,  it  stresses  using  experimental  and  theoretical 
research  on  the  mechanisms  of  raising  laser  damage  threshold 
values  to  draw  conclusions  and  make  inquiries. 

KEY  WORDS  Optical  coating  Laser  damage  Laser 
conditioning  Mechanism 
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1  INTRODUCTION 

Following  along  with  the  wide  spread  application  of  high 
power  lasers  in  fields  such  as  tactical  laser  weapons  as  well  as 
laser  thermonuclear  fusion,  and  so  on,  reserach  associated  with 
optically  coated  reflector  mirrors  possessing  high  laser  damage 
threshold  values  has  received  more  and  more  serious  attention 
from  people.  The  reason  is  that  laser  damage  threshold  values 
associated  with  optical  coatings  limit  laser  output  powers.  In 
the  last  more  than  ten  years,  a  new  type  of  technology  has  been 
developed — laser  conditioning  technology.  It  is  capable  of 
making  laser  damage  threshold  values  associated  with  certain 
highly  reflective  mirrors  go  up  two  to  three  fold.  Inside  China 
and  abroad,  a  number  of  hypotheses  have  been  put  forward  in  order 
to  explain  this  phenomenon.  What  is  regretable  is  that  there  is 
still  not  one  type  of  theory  capable  of  completely  explaining  it. 

2  CHARACTERISTICS  OF  LASER  CONDITIONING 

2 . 1  Laser  Conditioning  Technology 

As  far  as  laser  conditioning  is  concerned,  it  refers  to  a 
new  technology  which  is  capable  of  raising  2—3  fold  laser  damage 
threshold  values  after  taking  optical  coatings  or  large  pieces  of 
optical  materials  and  putting  them  through  laser  irradiation 
lower  than  laser  damage  threshold  values .  There  are  two  types  of 
laser  irradiation  methods  associated  with  laser  conditioning. 

One  is  using  multiple  iterations  of  laser  irradiation  of  the  same 
strength  (s-on-1).  Another  type  is  using  multiple  iterations  of 
laser  irradiation  with  optical  strengths  which  gradually  go  up 
over  time  (r-on-1).  Using  Hf02/Si02  and  Zr02/Si02  reflective 
film  systems  as  examples,  laser  damage  threshold  values  before 
laser  conditioning  are  7.4J/cm2  and  7.1J/cm2.  Going  through 
laser  conditioning,  they  rise  respectively  to  19.8J/cm2  and 
18.8J/cm2.  Moreover,  the  latter  type  of  laser  conditioning 
method  raises  laser  damage  threshold  values  somewhat  higher . 

2.2  Relationships  to  Laser  Parameters 

Laser  system  parameters  are  different.  Not  only  are  damage 
threshold  values  different,  but  damage  mechanisms  also  differ. 
Relationships  between  damage  threshold  values  and  wave  lengths 
have  roughly  the  relationship  •  x  is  a 

coefficient  which  changes  along  with  film  properties.  The  longer 
laser  wave  lengths  are,  the  higher  damage  threshold  values  then 
are.  The  relationship  between  damage  threshold  values  and  pulse 
widths  is  .  When  pulse  widths  are  long,  damage  is 

primarily  caused  by  absorption  of  focused  heat.  When  laser 
pulses  are  short,  dielectric  breakdown  created  by  electric  fields 
is  the  primary  factor  in  damage.  Besides  this,  there  are  also 
relationships  to  such  things  as  repetition  frequencies  of  lasers, 
modes,  and  so  on.  Due  to  the  fact  that  the  relationships  of 
laser  damage  threshold  values  and  widths  to  wave  lengths  are  the 
same  before  and  after  laser  conditioning,  as  a  result,  damage 
mechanisms  should  also  be  the  same  before  and  after  laser 
conditioning.  This  article  limits  itself  to  neodymium  glass 
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laser  spectral  lines  associated  with  irradiation  laser  light  of 
1.06  microns.  Pulse  widths  are  a  few  nanoseconds. 

2 . 3  Relationships  to  Methods  of  Coating 

Three  types  of  methods  are  used  to  carry  out  tests:  (1) 
commonly  used  electron  beam  evaporation;  (2)  ion  beam 
spattering  (IBS);  and,  (3)  plasma  plating  (PP).  Making  use  /298 
of  the  latter  two  types  of  methods ,  the  plated  coating  apertures 
are  few  (densities  are  high).  Mechanical  properties  and 
environmental  characteristics  are  good.  However,  experiments 
clearly  show  that  their  laser  damage  threshold  values  are  still 
lower  than  films  associated  with  electron  beam  evaporation.  As 
far  as  carrying  out  laser  conditioning  on  them  is  concerned,  PP 
film  laser  damage  threshold  values  are  not  raised  much. 

Moreover,  IBS  laser  conditioning  is  not  necessarily  effective. 
Generally  speaking,  laser  conditioning  is  only  effective  with 
regard  to  electron  beam  evaporation  coatings. 

2.4  Relationships  with  Thin  Film  Materials 

Thin  film  materials  associated  with  high  reflection  mirrors 
in  the  visible  and  near  infrared  often  use  pairs  of  Hf02/Si02, 
Ti02/Si02,  and  Zr02/Si02  to  form  reflection  film  systems. 

Because,  after  Ti02/Si02  plating,  there  is  a  requirement  for 
drying  in  oxygen  in  order  to  reach  complete  oxydation  so  as  to 
reduce  absorption,  this  is  not  appropriate  for  reflection  mirrors 
with  large  apertures.  Besides  that,  in  plating  processes, 
Ti02/Si02  reflection  mirrors  cannot  avoid  producing  intense 
absorption  of  low  valance  oxides  of  titanium.  As  a  result,  we 
only  select  Zr02/Si02  and  Hf02/Si02.  Making  use  of  laser 
conditioning  associated  with  pulse  widths  of  Ins,  laser  damage 
threshold  values  associated  with  the  latter  two  are  all  raised 
roughly  3  fold.  Replicability  of  laser  damage  threshold  values 
associated  with  Hf02/Si02  is  very  good.  Moreover,  levels  of 
increase  with  regard  to  various  types  of  pulses  are  the  same. 
Laser  damage  threshold  values  associated  with  Zr02/Si02  follow 
along  with  different  positions  on  reflector  surfaces  and  change. 
Laser  conditioning  results  are  sometimes  good  and  sometimes  bad. 
However,  making  use  of  Nomarski  microscopic  techniques  for 
observations,  the  flaw  densities  for  the  two  are  the  same. 

2.5  Relationships  with  Coating  System  Designs 

With  regard  to  three  types  of  samples — lainbda/4  Hf02,  rich 
Si02  samples  (includes  13  layers  of  lambda/4  Si02  between  lOnm 
Hf02),  and  rich  Hf02  samples  (includes  17  layers  of  lambda/4  Hf02 
between  lOnm  Si02) — after  laser  conditioning,  laser  damage 
threshold  values  all  go  up  approximately  2  fold.  Because  of 
this,  there  is  no  relationship  between  conditioning  results  and 
relative  thicknesses  of  film  systems.  However,  there  is  a  type 
of  exceptional  special  case.  That  is  nothing  else  than  that  when 
the  Si02  thickness  is  O,  that  is,  there  is  single  layer  Hf02, 
there  are  no  conditioning  effects. 
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3  INQUIRIES  INTO  LASER  CONDITIONING  MECHANISMS 
3 . 1  Early  Explanations  of  Laser  Conditioning 

The  first  type  of  hypothesis  is  "laser  cleanness", 
that  is,  laser  conditioning  takes  volatile  impurities  (in 
particular,  water  in  the  atmosphere)  and  gets  rid  of  them  from 
surfaces  and  the  interior  of  coating  materials.  These  impurities 
lead  to  laser  damage.  If  this  is  the  case,  damage  threshold 
value  increases  should  be  temporary.  We  took  Hf02  reflection 
mirrors  which  had  not  gone  through  laser  conditioning  and  put 
them  separately  into  environments  at  20 °C  with  relative 
humidities  of  approximately  50%  for  a  period  of  time.  Laser 
damage  threshold  values  of  1064nm  are  measured.  After  that,  use 
is  then  made  of  anhydrous  P205  and  measurements  are  made  again  in 
10  weeks.  Observations  are  done  of  the  shortening  of  wave 
lengths  after  the  removal  of  water.  It  was  discovered  that, 
comparing  laser  damage  threshold  values  before  and  after  the 
removal  of  water,  there  were  no  alterations  at  all.  However, 
reflectivity  drift  is  reversible.  The  explanation  of  this  is 
that  laser  damage  threshold  value  changes  are  permanent. 

Another  type  of  hypothesis  is  that  lasers  take 
materials  and  heat  them  right  up  until  their  crystalinities 
change.  In  order  to  check  this  hypothesis,  use  was  made  of  TEM 
and  X  ray  diffraction  in  order  to  determine  crystal  structures 
associated  with  single  and  multiple  layer  films  of  Hf02  and  Si02. 
Results  clearly  showed  that  Si02  in  single  and  multiple  layer 
films  of  1  micron  and  266  nm  was  no  different  at  all.  In  all 
cases,  it  was  in  amorphous  states.  In  contrast  to  that,  Hf02 
crystal  structures  were  related  to  thickness — smaller  than  20nm, 
states  were  indefinite.  As  growth  continued,  by  contrast,  they 
were  cylindrical  crystals.  They  were  primarily  monoclinic 
crystals.  Small  numbers  of  cubes  were  also  included.  However, 
changes  in  relative  thickness  certainly  did  not  influence 
conditioning  effects.  As  a  result,  there  was  no  relationship 
between  laser  conditioning  and  crystal  structures.  Making  use  of 
diffraction  meters  to  carry  out  scanning  with  regard  to  thin 
films  before  and  after  laser  conditioning,  no  differences  were 
discovered. 

3.2  Electron  Defect  Model 

Research  personnel  of  the  Lawrence  Livermore  National 
Laboratories  (LLNL)  put  forward  the  existence  [1]  of  a 
relationship  between  laser  conditioning  and  electron  defects 
inherent  in  films.  Because  electron  beam  evaporation  processes 
are  high  speed  quenching  from  high  temperature  ste^  to  low 
temperature  substrate,  as  a  result,  a  porous  fibn  is  produced 
with  defects  in  its  interior.  There  is  a  relationship  between 
nonequilibrium  states  and  laser  conditioning  characteristics. 
They  believe  that  laser  damage  associated  with  thin  films  goes 

through  the  four  steps  below:  ^  •  a.  j 

(1)  lasers  take  valence  carrying  electrons  from  restricted 

zones  and  excite  them  into  conducting  zones; 

(2)  through  photoelectric  field  or  free  current  carrying 
particle  absorption  effects,  free  current  carrying  particles  are 
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stimulated  to  achieve  high  energies; 

(3)  through  ionization  avalanche  or  photoelectric  effects 

(crystal  lattice  heating)/  excessive  amounts  of  energy  are 
transmited  into  crystal  lattices;  /299 

(4)  film  materials  are  heated  to  critical  damage 
temperatures — for  instance,  media  material  melting  points. 

The  occurence  of  laser  conditioning  phenomena  is  due  to  the 
elimination  of  the  sources  of  conducting  zone  electrons  in  the 
first  step,  that  is,  as  far  as  media  material  are  concerned, 
under  low  flux  laser  irradiation,  electrons  in  defects  are 
excited  into  conduction  zones.  However,  at  this  time, 
photoelectric  fields  are  still  not  adequate  to  give  rise  to 
damage  (steps  2-4).  Electrons  then  decay  into  relatively  deep 
energy  levels.  After  that,  it  is  not  easy  for  them  to  be 
stimulated  into  conduction  zones .  When  laser  conditioning  media 
make  use  of  laser  light  below  laser  damage  threshold  values  for 
irradiation,  the  number  of  electrons  it  is  possible  to  use  in 
order  to  transfer  to  conduction  zones  is  already  very  small.  As 
a  result,  the  amounts  of  net  energy  transmited  to  crystal 
lattices  is  also  small  and  does  not  go  so  far  as  to  cause  damage. 
In  this  way,  conditioning  effects  are  then  produced. 

Above  is  the  basic  thinking  associated  with  the  model  in 
question.  However,  LLNL  personnel  also  discovered  a  few 
experiments  [2]  that  do  not  agree  with  the  results  that  this 

model  predicts.  -oo 

Using  ND.YAG  lasers  with  wave  lengths  of  1064nm,  532nm, 

355nm,  and  pulse  widths  of  10ns,  respectively,  to  carry  out  laser 
conditioning  with  regard  to  1  micron  thick  Si02  and  Hf 02  single 
layer  and  multiple  layer  films,  it  was  discovered  that,  in  three 
wave  lengths,  Si02  single  layers  increased  roughly  2  fold  in  all 
cases.  However,  Hf02  single  layer  film  laser  damage  threshold 
values  still  showed  no  increase  at  all.  At  the  same  time,  LLNL 
personnel  determined  1.06  micron,  lambda/4  single  layer  Hf02 
Si02  films.  In  them,  Hf02  films  and  1  micron  thick  ones  are  the 
same.  There  are  still  no  conditioning  effects.  However,  the 
levels  of  increase  associated  with  Si02  laser  conditioning  laser 
damage  threshold  values  are  still  not  higher  than  30%. 

If,  with  respect  to  laser  damage  threshold  value  Dt^a.  ,  , 

X  =  approximately  1.0  for  Si02  single  films  and  x  =  approximately 
1.9  for  Hf02,  laser  damage  threshold  values  are  intensely 
influenced  by  wave  length.  In  another  area,  damage  forms  are 
different.  Si02  damage  is  a  few  small  points  (smaller  than  5 
microns),  and  they  do  not  get  any  bigger  with  continued  laser 
irradiation.  Yet  typical  Hf02  damage  has  diameters  or  0.1-1  mm. 
This  clearly  shows  that  the  methods  of  laser  coupling  to  Hf02  and 

coupling  to  Si02  are  different. 

With  regard  to  multiple  layer  Hf02/Si0  films,  after  1064nm 
laser  conditioning,  laser  damage  threshold  values  increase  2 
fold.  355nm  only  have  40%.  533nm  is  then  between  the  two. 

Laser  damage  threshold  values  at  1064nm  are  controled  by  Si02 . 
However,  at  very  short  wave  lengths,  by  contrast,  the  control  is 
by  Hf02.  Due  to  Hf02  absorption  fringes  associated  with  electron 
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beain  evaporation  at  around  375nin,  narrow  band  linear  absorption 
inside  Hf02  may  possibly  be  the  cause  controling  355nm  multiple 
film  layer  laser  damage  theshold  values . 

As  a  result,  it  is  certainly  not  a  single  electron  defect 
inside  Hf 02  or  Si02  which  is  then  capable  of  determining 
conditioning  effects.  In  that  case,  are  electron  defect  models 
right  in  the  final  analysis  or  not?  Up  to  now  it  has  still  not 
been  possible  to  reach  a  conclusion.  Perhaps  a  part  of  them  is 
correct. 

3.3  Recent  Progress 

Edwards  and  others  carried  out  analyses  [3]  with  regard  to 
Im  Zr02/Si02  Nova  laser  reflection  mirrors  plated  by  LLNL  (The 
surface  of  the  mirrors  in  question  is  unusually  smooth. 
Roughnesses  are  2.5-5nm).  They  discovered  that  the  production  of 
damage  is  begun  from  microscopic  defects  with  diameters  of 
approximately  30  microns.  These  defects  form  damage  of 
approximately  250  microns,  sometimes  covering  over  20%  of  the 
mirror  surface.  These  defects  lie  inside  thin  films.  In 
conjunction  with  this,  they  are  of  a  two  dimensional  "disc"  type. 
What  is  peculiar  is  that  the  damaged  film  surface  is  enveloped  by 
a  layer  of  "fog".  They  extracted  these  "fog"  granules  and 
discovered  that  they  were  precisely  Zr02  and  Si02  film  materials. 
The  explanation  for  this  lies  in  material  being  splattered  out  in 
the  damage  process.  In  conjunction  with  this,  it  piles  up  on  the 
film  surface.  Making  use  of  SEM  to  observe  these  granules, 
diameters  did  not  reach  1  micron.  There  was  both  surface 
smoothness  and  surface  roughness.  The  former  may  possibly  be 
formed  from  gas  states  condensing.  However,  the  latter  may 
possibly  be  solid  state  remnants  thrown  out. 

Edward's  explanation  is  that  it  is  due  to  the  defects  in 
Zr02  material  being  heated,  after  laser  irradiation,  to 
temperatures  sufficient  to  make  the  surrounding  Zr02  gasify,  that 
is,  in  excess  of  3000K.  In  reference  [4],  Chase  confirms  that 
the  cause  which  leads  to  this  process  is  linear  absorption.  Is 
this  type  of  so  called  defect  a  pure  substance  granule  or  not? 
Schildbach  made  use  of  laser  mass  analysis  systems  (LIMA),  and, 
in  conjunction  with  this,  discovered  no  impurities  at  all  [5]. 

In  another  area,  due  to  defect  densities  being  on  the  order  of 
lOOmm-2  (with  regard  to  the  Nova  reflection  mirror),  it  is  only 
when  a  vaccuum  chamber  is  severely  polluted  that  it  is  possible 
to  create  this  kind  of  outcome.  In  that  case,  may  they  not  be 
low  valence  oxide  materials — for  instance,  Si02?  Due  to  the  fact 
that  the  Nova  mirror  laser  damage  threshold  value  is  7.8J/cm2,  in 
this  way,  absorption  coefficients  associated  with  low  valence 
oxides  should  be  greater  than  lOOOcm-1  .  Unfortunately,  Si02 
absorption  coefficients  (1  micron)  still  do  not  reach  lOOcm-1  . 
What  defects  in  Si02  are  in  the  final  analysis  remains  a  puzzle 
even  now.  /300 

With  regard  to  Hf02  films,  LLNL  personnel,  when  studying 
electron  paramagnetic  defects,  discovered  that,  contained  in  rich 
Hf02  samples,  there  were  special  signals  which  were  not  in  other 
samples.  This  type  of  signal  clearly  shows  that  Hf02  may  contain 


32 


in  it  Fe3+  impurities  (Commercially  used  Hf02  film  materials 
contain  less  than  0.05%  Fe) .  This  cannot  be  eliminated  under 
laser  irradiation  below  damage  threshold  values,  leading  to  there 
being  no  conditioning  effects . 

Schildbach  discovered  that,  when  laser  irradiation  lower 
than  the  laser  damage  threshold  (16-18J/cm2)  is  used  on  Hf02/Si02 
reflective  films,  there  are  neutral  particles  of  Hf,  Si,  and 
oxygen  ( including  0  and  02 )  sent  out  [ 5 ] .  During  initial 
emission,  laser  optical  strength  is  10J/cm2.  Moreover,  the 
critical  optical  strength  required  for  conditioning  is  also  just 
10J/cm2.  As  a  result,  the  emission  of  neutral  particles  is  very 
possibly  produced  by  laser  conditioning.  Neutral  particle 
kinetic  energies  reach  as  high  as  a  few  score  eV.  High  energy 
neutral  particle  emission  corresponds  to  heat  induced 
f ractomission .  Even  if  there  is  no  obvious  damage,  this  may 
still  be  the  source  of  neutral  particles.  After  heating  of  thin 
film  defects,  heat  induced  cracking  is  produced,  leading  to  the 
elimination  of  absorbed  material  or  the  improvement  of  thermal 
stress  fields  (These  are  all  causes  giving  rise  to  damage). 

At  the  same  time,  Schildbach  analyzed  8  points  at  different 
locations  on  mirror  surfaces.  He  discovered  that  the 
conditioning  results  associated  with  them  were  certainly  not  the 
same.  With  regard  to  the  2  points  which  had  not  gone  through 
processing,  laser  damage  threshold  values  were  18“20J/cm2 — 
approaching  laser  damage  threshold  values  associated  with  other 
points  after  processing.  He  believed  that  the  cause  for  this  was 
due  to  these  two  points  not  having  any  defects. 

Staggs  and  others  used  atomic  force  microscopes  (AFM)  to 
observe  reflection  mirrors  after  plating.  They  discovered  that 
most  Si02  surfaces  contain  small  mounds  200nm  wide  and  lOnm  high 
[6].  Moreover,  Hf02/Si02  reflective  film  build  ups  are  primarily 
irregular  and  are  shallow  pits  or  cores.  Longitudinal  dimensions 
are  1-10  microns.  Transverse  dimensions  are  approximately  1/3  of 
longitudinal  dimensions.  Densities  are  on  the  order  of  lOmm-2. 
After  damage,  the  largest  changes  are  the  emission  of  cores.  A 
round  tunnel  remains .  Staggs  paid  attention  to  their  being 
cracking  on  the  edges  of  the  tunnels.  The  cores  are  very 
sensitive  to  damage.  When  cores  grow  on  thin  films,  due  to 
shadow  effects,  moreover,  the  peripheries  of  cores  are  made  to 
produce  gaps,  and  they  become  flexible.  Another  characteristic 
of  damage  is  the  enlargement  around  8  fold  of  small  mounds  in 
areas  without  defects  in  Si02  surfaces,  that  is,  their  growth  or 
agglomeration.  This  type  of  surface  agglomeration  temperature  is 
high  enough  to  reach  an  ability  to  make  components  transition 
within  time  periods  on  the  order  of  nanoseconds,  causing  the 
production  of  structures  with  large  volumes  and  small  surface 
areas.  Structures  clearly  show  that  emission  of  cores  occurs 
when  using  laser  irradiation  lower  than  laser  optical  strengths 
associated  with  making  small  mounds  agglomerate. 

Staggs  carried  out  irradiation  on  defect  areas.  During  the 
process  of  optical  strengths  increasing  from  7.8J/cm2  to 
54.6J/cm2,  AFM  was  used  to  analyze  the  appearance  of  defects. 
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There  were  still  no  obvious  changes  to  the  original  outlines. 
However,  after  increasing  to  19.8J/cin2,  most  lOnmxlOOnm  cores 
disappeared  from  surfaces.  After  increases  continued  to 
54.6J/cm2,  there  were  no  further  changes.  That  is  nothing  else 
than  to  say  that,  as  far  as  the  production  of  surface  smoothness 
phenomena  in  optical  strengths  emited  by  cores  is  concerned,  even 
though  optical  strengths  continue  to  increase — exceeding  laser 
damage  threshold  values — at  the  locations  associated  with  these 
cores,  still  no  damage  occured. 
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4.  CONCLUSIONS 


To  suiimarize  what  was  discussed  above,  relfection  mirror 
thin  film  damage  is  capable  of  definitely  not  resembling  that 
type  of  complete  cause  in  microscopic  nonstable  states  considered 
by  electron  defect  models,  and,  at  the  same  time,  is  also  not 
capable  of  being  explained  in  all  phenomena  from  a  macro 
viewpoint.  Combining  electron  defect  models,  we  believe  that  the 
damage  may  go  through  the  several  steps  below: 

(1)  Defects  Inside  Thin  Films — Under  Laser  Irradiation — 
Stimulate  the  Production  of  Large  Amounts  of  Free  Electrons 
Method  A:  Electron  Ionization  Avalanche 
If  free  electrons  which  have  broken  out  of  their  bonding 
already  exist  inside  thin  films,  after  laser  irradiation,  they  go 
through  electron  ionization  avalanche  and  form  plasma  bodies. 

The  avalanche  speed  is  represented  as 


dN, 

dx 


=  y{E)^N, 


In  the  equation,  y(E)  is  ionization  avalanche  coefficient,  Nc  is 
conductance  belt  electron  density,  and  E  is  field  strength. 

After  that,  energy  absorption  proceeded.  Assuming  that 

is  the  number  of  electrons  with  energies  between 
.  It  is  expressed  as 

»?(<J ,  T)  =  f/(Oexp((jXT) 
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Moreover,  electron  flux  through  isoenergetic  planes  is 

In  this,  v(<^).  and  U(0  are,  respectively,  effective  speeds 
and  coefficients  of  diffusion  associated  with  electron  flows. 
Then,  the  energy  distribution  equation  is 


^  =y(£)x»j(0 

dQ 


Attached  boundary  conditions  are 

,  t)=0 

J(0  ,t)=2xJ(£,,t) 

In  equations,  Eg  is  the  electrical  gap  of  materials. 

Method  B:  Going  Through  Valence  Electron  Transitions  Given 
Rise  to  by  Multiple  Photons 

At  this  time,  the  free  electron  densities  produced  by 


stimulated  valence  electron  transitions  are 


dx 


In  the  equation,  N  is  particle  density,  am  is  multiple  photon 
absorption,  Fm  is  incident  photon  flow  density. 

(2)  Free  Current  Carrying  Particle  Absorpton  Inside  Defects 
If  crystal  lattice  temperature  is  T.  On  the  basis  of  Drude 
models 


dT  N,  ehE^ 


dx  ni(l 

In  the  equation,  c  is  specific  heat,  p  is  crystal  density,  t  is 
photon  and  electron  collision  time  (approximately  on  the  order  of 
10-16  s),  E  is  electric  field  strengths  (effective  value),  m  is 
electron  effective  mass,  and  <o  is  laser  angular  frequency. 

Combining  the  energy  distribution  equations  above,  the 
carryS^out  of  integration  over  the  entire  energy  space  give 
laser  energy  absorption  as 


dT  _ 


m 


Tn  the  eauation.  Lac  is  the  average  free  travel  associated  with 
eLc^Jo^nd  pStton  collisions,  vs  is  vertical  wave  speed,  and  m* 

is  electron  belt  mass.  ^  ^  ^  ^ 

HeatincTf  Hicrh  Density  Defects  Produce 

Agglomeration.  In  Conjunction  with  This  Cracks  which  Are  Made  to 

Cover  Thin  Films  Lead  to  Damage.  i 

Laser  conditioning  is  then  the  production  of  ij. 

thin  films  under  low  flux  laser  irradiation,  and,  in  coniunction 
with  that,  the  causing  of  neutral  particle  emission.  The  J^®sults 
lead  to  laser  damage  cores  being  ejected.  Conditioning  res 
depend  only  on  defect  properties.  Moreover,  there  is  no 
relationship  to  such  factors  as  thin  film  density,  layer  number, 

^The^sianificance  of  laser  conditioning  research  lies  not 
only  in  being  able  to  make  our  advances  and  improvements  in  laser 
cSidlSonlS  technology,  which  is  not  ^ature  It  ts  also 
being  able  to  make  us  advance  our  understanding  of  the  ba 
nature  of  laser  damage. 
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